Kondo screening by the surface modes of a strong topological insulator 

E. Orignac 
Laboratoire de Physique de I' Ecole Normale Superieure de Lyon, 
CNRS UMR5672, 46 Allee d'ltalie, F-69364 Lyon Cedex 7, France 

S. Burdin 

Univ Bordeaux, LOMA, UMR 5798, F-33400 Faience, France and 

CNRS, LOMA, UMR 5798, F-33400 Faience, France 

(Dated: April 26, 2013) 

We consider a magnetic impurity deposited on the surface of a strong topological insulator and 
interacting with the surface modes by a Kondo exchange interaction. Taking into account the 
warping of the Fermi line of the surface modes, we derive a mapping to an effective one dimensional 
model and show that the impurity is fully screened by the surface electrons except when the Fermi 
level lies exactly at the Dirac point. Using an Abrikosov fermion mean-field theory, we calculate the 
shape of the electronic density Friedel oscillation resulting from the presence of the Kondo screening 

V-i cloud. We analyze quantitatively the observability of a six-fold symmetry in the Friedel oscillations 

£U for two prototype compounds: Bi2Se3 and Bi2Te3. 
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Topological insulators are a new class of materials, with an insulating bulk and a conducting surfaceP^. The 
existence of four topological invariants guarantees the stability of these surface states against perturbations.SHSl Fur- 
^_» thermore, the topological invariants permit to distinguish two types of topological insulators, the weak topological 

insulators and the strong topological insulators. In a strong topological insulator, the surface modes at low energy 
form an odd number of Dirac cones. This is to be contrasted with a strictly two-dimensional conductor such as 
graphene-^, where the Nielsen-Ninomiya theorenP only permits an even number of Dirac cones: the presence of the 
gapped bulk is crucial to the formation of an odd number of Dirac cones. Experimentally, the material Bi x Sbi_ x 
is a strong topological insulator with t hree Dirac conea^ while Bi 2 Se 3 , Bi 2 Te 3 and Sb 2 Te 3 , are strong topological 
insulators with a single Dirac conelMJJTJln t ne case f Bi 2 Te3, a significant hexagona l Ferm i line warping is present.^ 
More recent examples of topological insulator materials are TlBiSe^^l, strained HgTe^E^ Bi 2 Tei.6Si.J^PbBi 2 Tei^, 
i ^, Pb(Bii_ x Sba;) 2 TeJ^l, Bi 2 Se 2 Te, Bi 2 Te 2 Se, and GeBi 2 Te4p21 It has also recently been proposed on the basis of band 
structure calculations that the ternary compound LiAuSe^ and the cerium filled skutterudites^ CeOs4Asi 2 and 
CeOs4Asi 2 , could be topological insulators. A recent review of known topological insulators can be found in Ref. [STJ 
In all these systems, a strong spin-orbit coupling is present, and in the surface states the helicity, i.e., the sign of the 
spin projection of the quasiparticle spin on the quasi-momentum is fixed. An interesting theoretical question is then 
whether the Kondo effeclp2H2S] m the surface modes is affected by a fixed helicity. In particular, one would like to know 
whether a conventional Kondo effect takes place, or whether the fixed helicity constraint gives rise to unconventional 
Kondo fixed points. In the case of two-dimensional topological insulators^, the question was addressed^, and it 
was shown that a conventional Kondo effect would take place, leading to a suppression of the backscattering of the 
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edge states by magnetic impurities. In the case of a three-dimensional topological insulator, the Anderson model 
was considered both in the absence^®! anc i i n the presence^ of Fermi surface warping. Within a variational trial 
wavefunction method, it was found that the local moment would be fully quenched, but correlations would exist 
between the conduction electron spin and the spin of the local impurityPS Analogous results were obtained in a two- 
*\^ dimensional electron gas with Rashba and Dresselhaus spin-orbit couplings !^1 In the case of a pure Dirac spectr um ? a 
mapping on the one-dimensional Anderson model with a pseudogap in the hybridization function was obtained! 29 ^ 30 ! 
It was concluded that away from the Dirac point the Kondo effect would take place, and the impurity would be 
fully screened, while at the Dirac point, the local moment would decouple p2 The local density of states (LDOS), the 
local spin density of states (LSDOS) and the Friedel oscillations were also investigatecP^l within a U(l) slave-boson^ 
mean-field theory. In Ref l31l the Anderson model in a topological insulator with a Fermi surface warping was con- 
sidered using the numerical renormalization group and the beh avior of the LDOS was obtained. Experimentally, 
magnetic impurities such as manganese^, nickel, irorJ2£H23 cobaltP^EH and gadolinium 39 have been deposited on the 
surface of topological insulators. It was found that the surface states were remarkably insensitive to the presence of 
both magnetic and nonmagnetic impuritie d 34 * 39 l While the first result can be understood as a consequence of the 
suppression of backscattering, the second result is more surprising since magnetic impurities permit backscattering 
by flipping the electron spin. 

Besides this single impurity behavior, it has been suggested theoretically that in a Kondo lattice at electronic 



half-filling, the Kondo interaction could help the form ation of a topological insulator.ESHMl There are indeed recent 
experimental indications that the Kondo insulator 44 45 SmB 6 could be a topological insulator. This also lends support 
to the hypothesis that Kondo screening of magnetic impurities is compatible with the helical character o f the surface 
states. An important technique to probe conducting surfaces is Scanning Tunneling Microscopy (STM) SSIiZI This is 
particularly interesting in relation to the Kondo effect since STM measurements of the LDOS around a Kondo impurity 
located on the surface of a metal have already been performecP^SJ, anc j the influence of the Kondo resonance on the 
LDOS has been studied theoretically^!. Since the surface of a topological insulator is conducting, it can be probed by 
STM!3HE3The existing predictions for the LDOS caused by a Kondo impurity 3 -^ could thus be tested in that manner. 
Moreover, following the proposal of Ref. E2 integrating the measured LDOS would permit the measurement of the 
Friedel oscillations of the electron density induced by the Kondo screening cloud around the magnetic impurity. 

In the present paper, we want to further analyze the Kondo effect of a magnetic impurity at the surface of a strong 
topological insulator with warping. In the first part, we reduce the Kondo Hamiltonian to a one-dimensional model 
which can be treated by integrability techniques. We find that a conventional Kondo effect is obtained, with the 
impurity screened by the surface modes unless the Fermi level is right at the Dirac point, in which case, because of 
the vanishing density of states, the impurity remains unscreened for weak Kondo coupling. In the second part, we 
calculate within an Abrikosov fermion mean- field theory^! the Friedel oscillations 61 resulting from the existence of the 
Kondo screening cloud. For weak Fermi surface warping, we derive a perturbative expression of the density Friedel 
oscillations. In the third part we discuss the observability in STM measurements of the Friedel oscillations of electron 
density^!] around an impurity in the specific cases of two prototype compounds: Bi2Se3 and Bi2Te3. 

II. MAPPING TO A ONE-DIMENSIONAL MODEL 

The free electrons Hamiltonian of the surface modes of a strong topological insulator with warping reads: 

H =^2 clj-iv F i-{a a pxk) + -l(k x +ik y ) 3 + (k x -ik y ) 3 ]<j z - f i5 al3 ]c K i 3 = ^ 4,a W °/s( k ) c k./9 , (1) 

k,a„S k,a,/3 

where Ck, a annihilates a fermion of momentum k and spin a, z is the normal to surface of the topological insulator, 
a = (cr x , a y ,a z ) denotes the Pauli matrices, vp is the Fermi velocity, A is the war ping, a nd /i the chemical potential. 
In the case of Bi 2 Se 3 , one hasp! Vp = 3.55 eVA and A = 128 eVA 3 . For Bi 2 Te 3 J±ifei] one has v F = 2.55 eVA and 
A = 250 eVA. In the following we will study the effects induced by the presence of a magnetic quantum impurity. 
We take the position of the impurity as the origin of our coordinates so that the Kondo Hamiltonian describing the 
surface modes and their interaction with the impurity reads: 

H = H + H K = H Q + ^ J2 S ' c L^/3C k ^ , (2) 

where Jr- is the Kondo interaction, L the linear dimension of the (square) surface, S the impurity spin. 

A. Eigenstates of the free Hamiltonian 

In the present section, we briefly review the eigenstates of the free electrons Hamiltonian (TT1) . We will use a spinor 
notation to represent the Fermion annihilation and creation operators: 

*w = (^5) = E/^*w^). (3) 

\ / s—± 

where the spinors V'k( r ) are eigenstates of the first quantized Hamiltonian and {c s (k),<v(k')} = (27r) 2 (5 s , s '(5(k — k'). 
Introducing the polar coordinates k = (k, 4>) and r = (r, 9) the spinors have the explicit form: 



where \A S \ 2 + \B S \ 2 = 1 and 
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A s + iB s Xk 3 cos W + ivFk ._, 

s — = . (5) 



A s - iB s y /v 2 F k 2 + X 2 k (i cos 2 3e 



We note that A s and B s are periodic functions of 6 of period ^-. In the rotationally invariant case, i.e., without the 
warping term (A = 0), t he eigenstates of the Hamiltonian are rewritten as angular momentum eigenstates, and only 
the s-wave channeP2E2l is found to contribute to the Kondo interaction. In the case with warping (A ^ 0), there is 
only a discrete rotational symmetry, and instead of representations of 0(2) the eigenstates have to be expressed as 
representations of the discrete group Z^. 

«, W W4 E ^ £ "^ + ^(r), (6) 

* n=0,l,2 

with the restriction \<fi\ < tt/3 and £ = 0, 1, 2 labels the representation. With full rotational symmetry, £ would be the 
angular momentum. Because of the 3-fold symmetry of the warping term, the states having a difference of angular 
momentum equal to a multiple of 3 are hybridized together, and the angular momentum £ is defined only modulo 3. 
In that basis, the spinor ^ reads: 



*W= E /fC^MMMk) 



(7) 



Since ift%. ± 1 (0, 0) = 0, only the £ — 0, 2 modes interact with the magnetic impurity: 



, , d 2 k / Ac ,+ {k) ~ Bc a ^(k) . 

*(°) - ^ 3 / (2^)2 [ -*£c 2 , + (k) - a C2 ,_(k) I ' (8) 



and the Kondo Hamiltonian is rewritten as 
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// = E/(^( si? ( k )-^ c L( k ) c ^( k ) + ^ s -* t ( )^( ) . ( 9 ) 



(27T) 

with 

£(k) = ^v 2 F k 2 + \ 2 k 6 cos 2 (3<p) . (10) 

The Hamiltonian |9]) can be further reduced by turning the integration variable k to a system of curvilinear coordinates 
(E, k_l), where E is the energy of the eigenstate, and kj_ is the curvilinear coordinate along the line of constant energy. 
Introducing the new operators: 

r-c \ Q,s(k) 
ais(E, kj_) = . , (11) 

with anticommutation relations {ae, s (E, K±), a\ S (E', k'jJ} — (2tt) 2 S(E — E')8(kj_ — k'±), the free electrons part of the 
Hamiltonian Q can be rewritten: 

Ho = Z Zj d ^( sE -MjE,n ± )a e jE,n ± ) . (12) 

Introducing the density of states: 

»'W ^n£w (13) 



and the operators: 



T 7 v/2^Po(£)||V k £(k)|| 

i ^Vo(£)||V k £(k)|| 

a, + (£7) = ^ / ^ y ^K+( M. , (i 6 ) 

J V2* 2 A)(£)HV k £(k)|| 

a,,-(E) = ^/^^^ a2 -^^ , (17) 



with anticommutators {a a . s (E), a a , s i(E')} = (2n) 2 S(E — E')5 ss >S aa i we can define: 

a a (E) = {-i) l > 2 - a Y, &{sE)s 1 ' 2+a a a , s {\E\) , (18) 



s=± 



ip a {0) = / -= a a {E) , (19) 
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where denotes the Heaviside function. In terms of these operators, the Kondo Hamiltonian ffity reads: 

/+OO in 

j—^ Y( E - u)al{E)a a {E) + . . . , (20) 

where . . . stands for the modes of the free electrons Hamiltonian Hq that do not couple to the impurity. Away from 
the Dirac point, the density of states po(E) can be approximated by the density of states at the Fermi level pa(p), 
and the usual single channel Kondo problem is obtained. For a spin-1/2 impurity, The Kondo temperature can be 
obtained from the Bethe Ansatz solution as^ 

r, = De ™->M«), (21) 

where 7^ ~ 0.577 is the Euler-Mascheroni constant and D is a symmetric cutoff around the Fermi energy. The density 
of states po(E) can be expressed in terms of elliptic integrals 66 , so the full dependence of the Kondo temperature 
on the chemical potential is known up to the prefactor D. Since the Bethe Ansatz solution of the Kondo problerrPS 
requires a constant density of states, D represents as the energy scale away from the Fermi energy for which the 
density of states starts to deviate significantly from the density of states at the Fermi energy. For small warping A or 
not too far from the Dirac point, the density of states is a linear function of energy, and D ~ p, so the dependence 
of D on E is only a subdominant contribution. Close to the Dirac point, p, — > 0, the density of states po(p) ~ \p\. 
Because of such pseudogap, the Kondo t empe rature vanishes^. These results are in agreement with the ones derived 
in the framework of the Anderson modeP^ESl where Kondo screening was obtained only when the density of states at 



the Fermi level was non- vanishing. In the case of an impurity with spin S > 1/2 Eq. ( 20 1 would give the underscreened 
single channel Kondo fixed point. 

More generally, it can be established in any dimension d that with any free Hamiltonian of the form: 

Ho = Ycl a (p5 a p + e(k) • er Q/3 )c fei/ 3 , (22) 

a.fi 

having time reversal symmetry (i. e. e(— k) = — e(k)) only a conventional Kondo effect can be obtained. Indeed, if 
we write the partition function aipSl 

Z = Zo / Tre -JKj "drS(r)^i(0,r^ afl i, )3 (0,T)\ ^ (23) 

with Zq = Tr e~@ H °, and expand in powers of Jk, the series will depend on: 

" ( ) = ~p\:J^y (*n+/o 2 -Mk)ip e ' (24) 

but since e(k) is odd, introducing the density of states po(E) = J S(E — ||e(k)||)d' i k/(27r) d we can write: 

2/3^-^7 iv n + p, - E 

showing that the partition function is the same as the one of a system without spin-orbit coupling having the same 
density of states po a s the Hamiltonian (f22J) . As a result, a conventional Kondo effect is realized every time the density 
of states at the Fermi level is nonzero^. As a consequence, the dispersion and the Kondo self-energy are expected to 
remain spin-independent much beyond the mean-held approximation that we will consider in the following. 



III. FRIEDEL OSCILLATIONS AND ABRIKOSOV FERMIONS MEAN-FIELD THEORY 

A. Mean-field theory 

1. Abrikosov fermions and mean-field self-consistent relations 

We have seen in Sec ITT] that even in the presence of warping, a magnetic impurity on the surface of a topological 
insulator is always screened provided the density of states at the Fermi level is nonzero. In such conditions, a Kondo 
screening cloucF^ is formed around the impurity, and Friedel oscillations are formed. We consider in the following the 
case of a fully screened 5 = 1/2 impurity Since we are dealing with a conventional Kondo fixed point, we can use 
the Abrikosov fermion representation^! for the localized spin: 

S + = flh; 5- = /|/ t ; ff* = i(/f/t-/!/0, (26) 

with the constraint 

i = f\.h + f\h . (27) 

to rewrite the Kondo interaction as a local fermion-fermion interaction: 

H K = J K S ■ ttt( ) g *(0) = f f E tiWlMMO) - E €( °' 2 i(0) ) > (28) 

V a/3 a J 

where the two terms on the r.h.s. respectively correspond to the spin-flip and charge potential scattering processes 
of conduction electrons on the Kondo impurity Hereafter we will concentrate on the spin-flip interaction term and 
we will neglect the charge potential scattering one. Adding a Lagrange multiplier Hf(l — J2 a fafa) to the full Kondo 



Hamiltonian H + Hx to enforce the constraint (|27|) , we decouple ( 28 ) by a mean- field approximation: 



Hmf = H + A E(/1^(0) + H.c.) - n f E flfa, (29) 

a a 

where the effective hybridization A and the Lagrange multiplier jj,f satisfy the self- consistent relations : 

A=^E^l/a>, (30) 

a. 

i = E</«/«>' ( 31 ) 



where (• • •) denotes the thermal average computed with the mean-field effective Hamiltonian (29 1. We introduce the 
Fourier decomposition: 

^ Q (r) = ^E c k« elkr ' (32) 

k 

where L 2 is the surface of the system. To solve the mean field equations, we introduce the Green's functions 

G^(k,k',T) = -(T T c ka (r)c^(0)) I 
G^(k',r) = -(T T / a (r)c^(0)), 

G^(k,r) = -(T r c kQ (T)/t(0)), 

G f J p {r) = -(T T f a (r)f}(0)) . (33) 



Green's functions in ( 33 1 are expressed as 



Using the equations of motion from the Hamiltonian (291 and a Fourier decomposition in Matsubara frequencies, the 



G cc (k,k'>„) = a k ,*Go(k,u/ n ) + ^G (\c,w n )G^(w n )G (^',w n ) , (34) 

G /C (k>„) = ~G"(iv n )G (k,iv n ) , (35) 

G c/ (k,^„) = -^G (k, J ^)G // (^ n ) , (36) 

G*t(iv n ) = [iv n +n f - Z{iu n ))-\ (37) 



with the free electrons Green function 



and the self-energy: 



G Q {\L,iv n ) = {iv n -H\\<L))- 1 , 



£^„) = |A| 2 |^Go(k^ n ). 



(38) 



(39) 



Hereafter, we introduce the non-interacting electronic density of states, po(E) = j S(E — £k)d 2 k/(27r) 2 , w here ek 
denotes the electronic eigenenergies. Remarking that "H°(k) has time reversal symmetry, and invoking a similar 
analysis as the one leading to Eq. (25), we find the following spin-independent expression for the self-energy: 
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Po(E) 



z + p- E 
Introducing the real and imaginary parts of the self-energy, Y> a p(E + i0 + ) = 6 a ^[T,'(E) + iS"(E)], with 

1 
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de Po (e)P.V. r _ 
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Z"(E) = -Tr\A\ 2 p a (E + p), 
the self-consistency conditions ( 30 ) and pi] ) read: 

r °° dE 



|A| 2 = J K 

1 _ 

2 " 



_„ (F] (E + p f )X"(E) 

» tt F[ '[E + p f -^(EW + [^(EW 
dE £"(£) 

V" f( '[E + n f -Z'(E)] 2 + [Z"(E)]2 ' 



(40) 

(41) 
(42) 

(43) 
(44) 



2. Kondo temperature 



On general grounds, the Kondo temperature Tk indicates a crossover between the high temperature weakly coupled 
and the low temperature strongly coupled regimes. Indeed, the system at temperatures below Tk is characterized by a 
magnetic confinement of the spin of the edge electronic states which screens the Kondo impurity 33 . Tk has been shown 
to be the unique energy scale that characterizes the universal physical properties of single impurity Kondo models 
at low temperature T <C Tk- We will see later how this scaling property will become extremely useful for analyzing 
universal properties of the electronic density at very low temperature. But before, we derive here an expression of 
Tk- Within the mean-field approximation, the Kondo crossover turns to a transition at Tk which corresponds to 
a continuous vanishing of the / — c effective hybridization: A(Tr-) = 0. Invoking the self-consistent relations (43) 



and (44), we find (see appendix A 1 ) the following mean-field equation for Tk'- 
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and 
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P.V. 



+°° , n F (E) , 
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E 



E) 



(45) 



(46) 



which is equivalent to the Nagaoka-Suhl equation derived in Refs.|69]and[70l Assuming an even free electro ns den sity 
of state, po(—E) = po(E), a general explicit expression of Tk was derived in the weak Kondo coupling limfpM^] 



T, 



ith 



A" 



2e 7B 
F K = -^(D + p)exp 



F K e 



-1/JkPoM 



dE p (p + E) - po{p) 



E 



Mm) 



(47) 



(48) 



where je — 0.577 is the Euler-Mascheroni constant, and D denotes the half-bandwidth of po- The mean-field 



expression (47) provides the usual non-analytic exponential term characterizing the Jk~ dependence of the Kondo 
temperature at small Kondo coupling JkPo(p) "C 1. The same exponential dependence also emerges from the 



mapping to the Bethe Ansatz solvable model (see Eq. 21). The prefactor Fk does not depend on Jk and is thus a 



pure characteristics of the non-interacting electronic system. Whatever the specific chemistry of the magnetic impurity 
is and whatever the microscopic details of its coupling with conduction electrons, the prefactor F K depends only on 
the energy structure of the non-interacting conduction electrons. 

We find that Fk ~ D if one assumes a constant density of states and a chemical potential close to the middle of the 
electronic energy band. The situation may become quantitatively different when the conduction electrons describe 
the surface modes of a topological insulator. Indeed, using the density of states p{E) = \E\/(2nv\) that characterizes 
a surface mode without warping, we find for ji > 0: 



T K 



2e 7E " 



/'" 



D + fi 



exp 



D 



2-KVp 



fi J K p 



(49) 



This small Jk asymptotic expression becomes singular at the Dirac point [i — 0, where the density of states vanishes 
linearly At the Dirac point, one has to start from the equation (46) for the Kondo temperature, which simplifies to: 



1 

~J~K 



where D is a bandwidth cutoff. For Jk < J K 
Ref. E3 For J K > J K , we find that: 



D 



T K 

TTVJr 
2nv% 



dE , / E \ 

7t tann 

2irv 2 F \2TkJ 
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cosh | — I 
V.21W 



(50) 

(51) 



, the Eq. ( 50 1 has no solution, in agreement with the results of 



T K 



D 

2 In 2 



1- 



'K 

Jk 



(52) 



indicating that regular Kondo effect can still be realized at the Dirac point, but with a Kondo temperature vanishing 
linearly with the Kondo interaction when Jk — ► Jk j hi agreement with the prediction of Ref. [67] for a linear density 
of states. 



3. T-matrix and local electronic density 



The local electronic density is defined as: 

P (r) = -ilm(Tr[G cc (r,r,z0 + )])=p ( M ) + <5p(r) . (53) 

From Eqs. (34) and (37) we derive the expression of the T-matrix^, which is defined by the relation : 

G(r,v',iv n ) =G {r-r',ip n )+G (r,ii> n )T(iu n )Go(-r',w n ) . (54) 

We find: 

T{iv n ) = |A| 2 G"(u/ n ) = |A| 2 (^ n +nf- S(z^))- 1 . (55) 

Invoking this expression of the T-matrix and considering that the Kondo self-energy H(iv n ) is diagonal and symmetric 
in spin components (see Eq. (40)), the local electronic density is then given by : 



| A| 2 ^ Tr [G (r, w n )G (-T, w n )] 



P 



iv n + (i f - £(«/„) 



(56) 



Assuming that the energy scale which characterizes the Kondo resonance is much smaller than the effective band- 
width of the non-interacting electrons, the T-matrix can be approximated as: 



T(iv n ) 



IAI 



i{v n +rsign(i/„)) 



(57) 



where T = 7r|A| 2 /?o(^). This expression of the T-matrix is expected to be valid beyond the mean-field approach that 
we are following here, since it is equivalent up to a Fourier transformation to the definitions Eqs. (5.38)-(5.39) on 
p. 112 in Ref. [Ml 



The expression of T in terms of the microscopic parameters of the Hamiltonian can be obtained from Eq. (A12| 
For zero temperature, we have (see App. [Alfor a derivation) : 



T(T = 0) = (D + fx) exp 



JkPo(p) 



dE p (p + E) - po(n) 



-D-fj, 



E 



Po{p) 



(58) 



Comparing this expression with expression (47) for the Kondo temperature, we find the very general (i.e, coupling 
and band-structure independent) relation 



2e 7E 

t k = — r(T = o) 



(59) 



which connects universally the crossover temperature Tk to the T = resonance width. Invoking the Wilson ratio 
R = x(T = 0)/7 = 2, and using the small coupling asymptotic expression 7 ps — ^Im (G^(z0 + )) ps 1/nT, we find the 
following universal Wilson number within the mean-field approximation: 



2e TB 
w = T KX (T = 0) = — 3- « 0, 36 



(60) 



This result is in relatively good agreement with Wilson's numerical result io = 0,41 (see Eq. (6.75) in Ref l3"3"|). 
Therefore, the following results will be derived within the mean-field method, but we expect them to be qualitatively 
and quantitatively valid beyond this approximation. 



B. Friedel oscillations without warping 



In this section, we study the electronic local density p(r) in the absence of warping, i.e., for A = 0. We start from 
expression ( 53 ) and we rewrite the sum as a contour integral: 



6 P (r) = -J c ^n F (z)Tr 



IAI 2 
G„(r,*)-i-i-Go(-r,«) 
z + 11 



c> 



dz 
2hr 



n F (z)Ti 



IAI 2 
G (r,z)^\G (-r 1 z) 
z — 11 



(61) 



where the contours C\ and G2, depicted by figure[l] encircle the poles of the function n F {z) = (e l3z + 1) _1 of positive 
(resp. negative) imaginary part. The contours can be deformed respectively into C[ and C 2 allowing to rewrite the 
integral as: 



— n F {E)lm 



Tr(G (r, E + i0+)G (-r, E + i0+)) 



E + iT 



(62) 



where E covers the real axis and Gg(r,E + i0 + ) will be obtained in the following from analytic continuation. In- 
deed, in the absence of warping, an analytic expression of the Green's function is available from the non-interacting 
Hamiltonian (nj). After Fourier transform of [iv n — %°(k)) -1 , we find : 



r< f • \ iv n+ P 

G {r,u> n ) = — — %- 
2irvt 



,. , p n \-iy sign(i/„) 
K r 

v F 



~, r . , Wn\ -ifJt sign(i/„) 
sign(zy n )z • (a x r)^ ( ^ — y -^r 



(63) 



and 



Tr[Go(r,iv n )G (-r,iv n )] 



(iv n + p) 2 
2tt 2 4 



j. , i^nj ~ip sign(^„) \ f \v n \-ip sign(^„) x " 



V F 



Vf 



(64) 



where K and K\ are respectively modified Bessel functions the second kind of order zero and onep^ and r = |r|. 
Hereafter we will replace this expression by its analytic continuation in the upper half complex plane, iv n — > E + i0 + , 
with u n > 0. Using Eq. (9.6.4) in Ref. UJ2 we have: 



K„ 



0+ -i{E + p) 



vf 



2 



j f(E + p + iQ+) 

VF 



(65) 



where Hi, denotes the Hankel function.^ We thus find the expression: 



Mr) 



871-4 



j dE (E + pf n F (E) Im 



(J^W) 2 (Mr) + (Ff 
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(66) 



For long distances, r 3> Rp = Vf/h, we can use the approximation from Ref. [75J 
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giving: 
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where we have used the relation |A| 2 = r/7rpo(/z) = 2Yv F / p. 
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FIG. 1. The integration contours C\ and Ci used in the calculation of the sum (53 1. These contours can be deformed into C| 
andCa. 



At zero temperature, with Eq.(68) we find the following expression for the local density: 
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where E\(u) is the exponential integraP^. For long distance r ^> vp/Y , since e u E\(u) ~ 1/u, we obtain 
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(69) 



(70) 



Remarkably, the amplitude of the Friedel oscillations become independent of T, i.e., independent of Tk, for distances 
longer than the Kondo lengthscale Rk = Vf/T. This can be understood by noting that for distances larger than Rk, 
the impurity appears as a potential scatterer at unit arity!^ For short distances r <C vf/(2/j,), Rk-, we find 



5p(r) 



T , [2Ye 1 r 
In 

v F 



2, /T ,2 



27r 2 /ir 



(71) 



so that the amplitude of oscillations inside the Kondo cloud depends explicitly on T. This indicates that the Friedel 
oscillations at distances shorter than the Kondo lengthscale Rk reflect the internal structure of the Kondo screening 
cloudP 



At finite temperature, Eq.(68) yields: 

r e"^ 



/ 1 BY 3 BY -ah 



(2pr\ 



(72) 



where 2-F1 is the Gauss hypergeometric function^. For long distances, and T -C Tk the Friedel oscillations decay 
exponentially over the thermal length R-p = vp/T, 
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(73) 
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FIG. 2. (Color online) Variation of electronic density computed at finite temperature using the exact expression Eq. ( 72 1 [red 
solid line], the long distance coherent approximation Eq. (|70| [blue doted line] which is appropriate for Rk < r < Rt, and the 
very long distance regime r > Rt approximated by Eq. ( |73[ ) [green dashed line], a) Overview along few Friedel oscillations, b) 
Focus at relatively short distance. Here, we chose Rt = 50Rk and irRp — 5Rk- 
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FIG. 3. (Color online) Pl ot o f the envelope of Friedel oscillations, <5p(r)/cos(2/ir/wF), compu ted at finite temperature using 
the exact expression Eq. ( 72 1 [red solid line] , the long distance coherent approximation Eq. ( 70 1 [blue dotted line] which is 
appropriate for Rk < T < Rt, and the very long distance regime r > Rt approximated by Eq. (731 [green dashed line]. Here, 
we chose Rt = 50Rk and ttR f = 5Rk- 



in agreement with Ref. I3TJ1 The behavior is represented on Fig. [2] which depicts the Friedel oscillations, and on Fig. [3] 
which depicts the envelope of these oscillations. It appears clearly that the Friedel oscillations have a period irRp, and 
an envelope decreasing like 1/r 3 when Rk < r < Rt as approximated by Eq. (70 1. Their amplitude is exponentially 
reduced at longer distances r > Rt where Eq. ( 73 ) provides a good approximation. We see from the previous 
analysis that there are three relevant lengthscales as depicted schematically by figure |4] The first one, Rf — vf/^, is 



proportional to the Fermi wavelength. For lengthscales much larger than Rp, the simplification (67) is justified, and 
Rp then simply gives the pseudo-period of the Friedel oscillations. The second lengthscale is the Kondo screening 
length Rk, which is the size of the Kondo cloud; due to the temperature dependence of the Kondo resonance width, 
r(T), we have Rk{T ^ Tk) — oo and Rk(T <§C Tr-) oc vp/Tk- Above that lengthscale, which requires at least 
T < Tk, the Friedel oscillations become identical to those of a resonant non-magnetic impurity. The third important 
lengthscale is the thermal length Rt- Beyond that length, the Friedel oscillations decay exponentially, while below Rt 
they are unchanged from the zero temperature case. In order to observe the Friedel oscillations characteristic of the 
Kondo screening cloud, we must have T ^C Tk so that the Kondo screening length is much shorter than the thermal 
length. In a renormalization group picture, the temperature is the natural infrared cutoff for the rcnormalization 
group flow, and the constraint T <C Tk is simply a requirement that the strong coupling scale is reached before the 
thermal cutoff. 
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FIG. 4. Schematic description of the Friedel oscillations that can be observed in the density variations around the impurity 



when T < Tk- For the rotationally invariant case A : 
is slower (like 1/r 2 ) in some directions. 



the envelope decays like 1/r . For the warping term A/0 the decay 



C. Friedel oscillations with warping 



We now turn to the electronic local density p(r) in the presence of warping. For A / we don't have anymore an 
expression in closed form of the Green's function. Instead, we use the approximation (|B13l derived in App|B] : 



with 



We obtain: 
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(77) 



Since this expression is only approximate, we cannot use contour integration techniques to obtain the sum ( 53 1 . Indeed, 



attempting to use contour integration for T > yields a divergent integral. Nevertheless, for zero temperature, we 
can still change the sum ( 53 ) into an integral as the exponential decay of the modified Bessel function ensures the 



convergence of the sums. This leads to the zero temperature result valid for r ^> Rk — vp/F 



Sp(r) = 5p(r)\ x =o 



^^ye )cos ( 2 pr\ 

v F J 



(78) 



For positive temperature, we have to compute the sum (53) numerically. This will be done in the next section using 



realistic values for the model parameters. Here we rather discuss general new features that emerge from this warping 
term. First, we expect finite temperature corrections to be relevant only in the crossover temperature regime around 
Tk- Indeed, similarly to what we found without warping, the thermal length Rt = vf/T provides a cut-off below 
which Friedel oscillations are identical to the ones predicted for T = 0, and above which they are muffled by thermal 
fluctuations. 

Furthermore, one important feature here in this expression is the 1/r 2 decay of the envelope: this decay is identical 
to the one of a two-dimensional normal metal, and it dominates the 1/r 3 contribution from the non- warping. Because 
of the cos 2 39 factor, the Friedel oscillations in the directions 9 — f (2n + 1) have the contribution from warping 
switched off and contain only the 1/r 3 contribution, whilst in other directions the warping contribution is observable 
and dominates on the longer lengthscales due to its slower 1/r 2 decay. 

Also, Rk = Vf/T characterizing the size of the Kondo screening cloud, we expect the Friedel oscillations to be 
observable only at distance larger than this Kondo length. Within a renormalization group picture, the density 
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n s w ng/10 



n s = 1,6.10 3 /znT 
/i = 50 meV 
Rf = 7, nm 
Rw = 130 /im 
T w = 3.10" ' 
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n s w 10n g 



Tk = 10 K 

i?K = 410 nm 



ns — 1, 6.10 4 um 
|U = 160 meV 
R F = 2,2 nm 
i?w = 410 nm 



n s = l,6.10 t> /im" 
u = 500 meV 
Rf = 7, A 
fliy = 1,3 nm 



Rk 



= 100 K 
= 41 nm 



4.10 3 urn -2 



ns 
P> = 

TV = 3.10~ 4 K 



80 meV 
4, 4 nm 
= 13 /im 



ns = 

l« = 

Rf 

Rw 



= 4.10 4 urn" 2 
250 meV 
= 1,4 nm 
= 41 nm 



ns - 



= 4.10 5 urn" 2 
800 meV 

= 4,4 A 
= 1,3 A 



Tif = 1000 K 
Rk = 4, 1 nm 



ns — 10 /Ltm 
u = 130 meV 
Rf = 2,8 nm 
Rw = 1,3 /im 
TW = 3.10" 3 K 



n s = 10 5 /im"' 
u = 400 meV 
Rf = 8,8 A 
_Rvk = 4, 1 nm 



ns = 10 /xm 
u= l,3eV 
i?F = 2,8 A 
i?w = 0, 13 A 



TABLE I. Model parameters for Bi 2 Se 3 , with v F = 3.55 eVA and A = 128 eVA 3 



oscillations are thus supposed to be measured at a sufficiently large distance from the Kondo impurity such that the 
system is correctly described by its strong coupling fixed point, i.e., the Kondo spin is totally screened. Nevertheless, 
comparing Eq. ( 70 ) and Eq. ( 78 ) we find that the warping correction becomes relevant only for distances larger 

— . . — v 7 

than a new characteristic length, Rw = A2 F 5 ■ Introducing the density of surface states ns and invoking the density 



of states po(E) — \E\/(2ttv f ) in the vicinity of the Dirac point, we have ns = p 2 /(Anvp) = l/(4irR F ), so that 

/ 2 \ 2/5 

Rw = (wfM) 2 (47™s)~ 5 / 2 . A crossover density emerges, n* s = -^ f A2 "^ j , that distinguishes two different cases: 

for ns > n%, we find Rw < Rk an d the Fricdcl oscillations which are observed for Rk < r < Rt are characterized by 
the warping term with switch on and off directions and an 1/r 3 envelope. But, closer to the Dirac point, i.e., for density 
ns < ng 7 Friedel oscillations are characterized by two regimes: at shortest distances Rk < r < Rw the isotropic term 
with 1/r 2 envelope dominates, whilst the warping correction dominates at larger distances Rw < r < Rt- Also, a 
new temperature scale emerges in the lowest density case ns < n^: the warping temperature Tyy = vp/Rw < Tk 
indicating the crossover temperature below which warping effects appear. 



IV. DISCUSSION 



A. Experimental observability of the density oscillations 



Here we analyze the experimental observability of the density oscillation effects, with or without warping effects, 
that were discussed on general grounds in the previous section. The idea is to compute the density variation Sp(r) 
around a Kondo impurity using realistic values of parameters that correspond to topologically insulating compounds 
for which surface states have been observed or predicted. We consider more specifically two compounds: Bi2Se3 for 
which Ref.foUl gives v F = 3.55 eVA and A = 128 eVA 3 , and Bi 2 Te 3 with values v F = 2.55 eVA and A = 250 eVA 3 given 
by Refs. ITT1 and [Ml We are still left with two tunable parameters: the Kondo temperature Tk = vf/Rk, and the 
density of surface states ns = p 2 /(Anvp) = l/(4irR F ). First, we remark that the Kondo temperature depends on 
various microscopic parameters including the chemistry of the magnetic impurity and the density ng. Furthermore 
the well know exponential Jr— dependence of Tk (see section III A 2 1 makes this temperature scale very sensitive 
to variations of these microscopic parameters. Therefore, refereeing from the orders of magnitude that are usually 
measured in Kondo compounds we consider here three different characteristic values: Tk — 1000 K (big), Tk = 100 K 
(medium), and Tk = 10 K (relatively small). For the density of surface states ns, we consider three values for each 

/ v 2 \2/5 

compound: the crossover density n* s = -^ I A2 "g j , a smaller density ns = rig/10, and a larger density ns = 10n s . 
The sets of relevant parameters that we consider are summarized in table [I] for Bi 2 Se 3 , and in table [TT] for Bi 2 Te 3 . 

The results are represented on Figs. [5] and [6] for Bi 2 Se3 and Bi 2 Te3 respectively. For these plots, we fixed arbitrarily 
Tk = 100 K and we chose realistic relevant values of chemical potential p, that can be controlled experimentally by 
doping with Sn 11 or Mg 63 . In Bi 2 Se3 compounds, experimental values of p indicated in Ref. 1 63 1 are tuned from 350 meV 
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TABLE II. Model parameters for Bi 2 Te 3 , with v F = 2.55 eVA and A = 250 eVA 3 



down to eV. Therefore, the four values that we considered for the plots of Fig.[5]were chosen invoking table[T]in order 
to illustrate the observability of the various cases: with dominant warping term (/i — 350 meV> /!*), with similar 
warping and isotropic terms (p = fi* = 250 meV), and with negligible warping term (/i = 130 meV, and 50 meV). 
Figure [5] clearly shows the Friedel oscillations with six- fold rotation symmetry when \x > fi*, or with full rotation 
symmetry when fi < fx* . The choice of Tk = 100 K for this plot is arbitrary and experimental values of the Kondo 
temperature can be significantly different. Furthermore, we are aware that doping, i.e., varying fi, strongly affects 
the value of Tk which may continuously vanish at the Dirac point as we discussed in section [ill A 2| Nevertheless, we 
expect that the Friedel oscillations will qualitatively not really depend on Tk- This is a consequence of the universality 
of the strong Kondo coupling effective regime that is realized below Tk within the renormalization group picture: 
since Friedel oscillations appear above the Kondo screening size Rk their shape is qualitatively universal (but still 
depends on the warping length Ryy and Fermi pseudo period irRp). Also, according to the values given in table [11 
the crossover value for the chemical potential varies very smoothly from (i* — 160 meV to fi* = 400 meV when Tk 
changes from 10 K to 1000 K. This suggests that the results which are illustrated by figure [5] for Tk = 100 K can be 
extended to any other values of Tk- Of course, the characteristic unit length which is used for the plots, Rk, would 
have to be rescaled accordingly. Experimentally, one of the main difficulties for observing Friedel oscillations with or 
without six- fold symmetry is the requirement of cooling the temperature sufficiently lower than Tk, but the orders 
of magnitudes that are considered here correspond to values that are realistic both physically (Tk is imposed by the 
chemistry) and technologically (T is limited by cryogenic technics). 

We made a similar analysis for Bi2Te3 compounds. In this case, Ref. [11] indicates experimental values of \x between 
350 meV and 120 meV. We thus plotted these two extreme cases, together with the intermediate value /i = 250 meV. 
Here, we chose a Kondo temperature Tk = 100 K, which corresponds to a crossover value fi* = 130 meV as indicated 
in table In] In this case, the six-fold symmetry resulting from the warping is thus observable for /i = 350 meV and 
250 meV, and the full rotation symmetry is recovered for // = 120 meV. The six-fold symmetry may remain for that 
value of chemical potential if the Kondo temperature is lowered. Indeed, table [TT] indicates fi* = 82 meV for Bi2Te3 
compounds with Tk = 10 K. 

Here, we have restricted our analysis to the observation of Friedel oscillations within the fluctuation of the local 
density of states, Sp(r). This physical quantity can be measured experimentally using Scanning Tunneling Microscopy 
(STM). Local Density of States^ (LDOS) measurements have already been performed by STM on Kondo impurities 
at the surface of metalsPSHSHThe measurement of the Friedel oscillations would require the integration of the measured 
local density of states over a range of energy^. Beside the issue of cooling the temperature sufficiently lower than 
Tk, other technical limitations have to be considered in order to observe the predicted Friedel oscillations using STM: 
First, a voltage bias is applied locally between the tip of the STM and the surface of the sample. The resulting STM 
current which is measured may invoke out of equilibrium effects that have not been analyzed here. We expect our 
predictions to be valid for STM experiments with bias voltages invoking energies that are lower to both Tk and /i. 
Higher values of bias voltage may have non universal effects on the Kondo screening leading to a distortion of the 
Friedel oscillations. 
A second limitation is the STM resolution in both lateral and depth directions. More precisely, we may expect an 
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FIG. 5. Intensity plot of the Friedel oscillations around a screened magnetic impurity at the edge of Bi2Se3. Model parameters 
v F = 3.55 eVA and A = 128 eVA 3 . We chose T K = 100A' and T = 10K. From top left to bottom right, /x = 350 meV, 
250 meV, 130 meV, 50 meV. Unit length: Kk- Short distances, r < 0.2Rk, are not represented. Intensity in arbitrary units is 
represented by the darkness of the plots. 



experimental STM measurement of the Friedel oscillations to be realized by moving the STM tip at the surface of 
the system along two orthogonal directions. The most natural resulting STM signal will thus be discretized on a grid 
having a square lattice symmetry and an elementary step of length Rstm — 0.2nm. Since the period of the Friedel 
oscillations is ttRf, the measured STM signal may exhibit a Moire patterrPH resulting from the interference between 
the two periods, Rstm and ttRf- Considering the values of Rp which are given in tables [T| and \U\ and assuming 
Rstm is of the order of one or few A, Moire patterns might occur for values of chemical potential relatively higher 
than il* . In such cases, the measured STM images would only have the two- fold symmetry common to both the 
square and the hexagonal symmetry groups. 

Comparing qualitatively the plots of figures [5] and [61 we find that warping effects and their related six- fold symmetry 
are more observable at the surface of Bi2Te3 rather than Bi2Se3. This is due of course to a larger value of the warping 
constant A, but this also results from a smaller value of the Fermi velocity vp, which gives a smaller value of crossover 
potential /i*. 
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FIG. 6. Intensity plot of the Friedel oscillations around a screened magnetic impurity at the edge of Bi2Te3. Model parameters 
v F = 2.55 eVA and A = 250 eVA 3 . We chose T K = lOO-ftT and T = IQK . From left to right, /j, = 340 meV, 250 meV, 120 meV. 
Unit length: Rk- Short distances, r < 0.2Rk, are not represented. Intensity in arbitrary units is represented by the darkness 
of the plots. 



B. Conclusion 



We have shown that a magnetic impurity on the surface of a strong topological insulator will be fully screened by 
the surface modes unless the Fermi energy is exactly at the Dirac point. The result depends only on the time reversal 
invariance of the effective Hamiltonian of the surface modes and is valid in particular in the presence of warping of 
the Fermi surface. We have shown that Friedel oscillations are formed around the impurity and we have calculated 
the shape of these oscillations both without warping and with a weak warping that can be treated perturbatively. 
With warping, the symmetry of the Friedel oscillation pattern is broken from full rotational symmetry to a six-fold 
symmetry. In both cases, the pseudo period of the oscillations, nRp = irvp/fi, is half the Fermi wave length of the 
surface modes, the short distance cut-off Rk oc vp /Tk is determined by the Kondo temperature Tk , and the long 
distance cut-off Rp = vp/T results from thermal fluctuations. With warping, the amplitude of the fully rotationally 
symmetric part decreases as 1/r 3 , whilst the 6-fold symmetry term has an envelope decreasing more slowly, a 1/r 2 . 
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As a consequence, a new length scale Rw emerges, above which Friedel oscillations with six-fold symmetry may be 
observed. The crossover condition Rk ~ Rw defines a chemical potential /i* associated to a doping rig above which 
the Friedel oscillations are characterized by the six-fold symmetry even at shortest distances. Considering realistic 
values for the model parameters, we analyzed the observability and the symmetry of Friedel oscillations in the vicinity 
of magnetic impurities deposited at the surface of two compounds, Bi2Se3 and Bi2Te3. We identified large range of 
parameters where the crossover between the 6-fold and the fully rotational symmetries may be observed. We propose 
to use STM as an experimental probe for the variation of the local density of states. 

Various questions remain to be addressed. First, it would be interesting to investigate the unstable fixed point 
that separates the regime with Kondo screening from the regime of decoupled impurities in the case of a system at 
half-filling. Second, an exact calculation of the Friedel oscillations could be performed using form factor expansion 
mcthodd 75 * 76 ! at zero temperature. 

The single impurity model that we considered here could also be generalized to two or several impurities. In this 
case, the local Kondo screening will compete with the Ruderman-Kittel-Kasuya-YosidcP^HISl (RKKY) inter-impurity 
screening as discussed by DoniacfP2 in a general context. Here, we may expect the 6-fold symmetry of the Friedel 
oscillations to have signatures on the symmetry of the RKKY interaction. 
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Appendix A: Details of the mean-field approximation 



We start from the mean- field self-consistent relations ( 43 ) 
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2 J_ x n [£ + /i/ -£'(£)]* + [E"(25)] S 



with the expression (42) of the imaginary part of the self energy: 

X"(E) = -tt\A\ 2 Po (u + E). (A3) 

Introducing the Hilbert transform of the density of states: 

*(E) = J~ depo(e)P.V. (^) . (A4) 



the expression (41 ) of the real part of the self-energy can be written as: 

£'(£)= 7r|A| 2 *(/i + £) . (A5) 

1. Derivation of the Kondo temperature 



When the temperature T goes to Tjc, the hybrid izati on parameter A vanishes. In such limit, the Eq. ( A2 1 reduces 
to n F (-/i F ) = 1/2, yielding Hf{T K ) = 0. The Eq. (|aT|) can then be cast into the Nagaoka-Suhl forirPOlU 
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In the general /i > case, we can rewrite Eq. ( A6 ) : 
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where D is a bandwidth cutoff such that p(E) = when \E\ > D. taking the limit of low Kondo temperature 
Tk <C D,p, we obtain the following approximate equation: 



1 (2ei»{D + n) 



D-fi 



dE 



[p (p + E) - p (p)) , 



(A8) 



improving the prefactor in the expression of the Kondo temperature by taking into account the variation of the density 
of states with the energy. 

2. Derivation of the resonance width 

We have seen that both T and pj vanish at and above the Kondo temperature. Hereafter, we will consider the 
limit of small Kondo coupling, and we will thus assume that T and pf remain small compared to the non-interacting 
electron characteristic energy scales, even below the Kondo temperature where these quantities are not vanishing any 
more. If we take first the equation (A2|, and make the approximations Y/{E) = 0, £"(1?) = E"(0) = — T, we obtain: 
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ir (E + p f ) 2 +T 2 
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(A9) 



so we have to take also pt = 0. Inserting our approximations in the denominator of ( Al ), we obtain a second equation: 

(A10) 
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We introduce the bandwidth cutoff D such that po(E) = if \E\ > D. We can then rewrite Eq. (|A10) in the form: 
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(All) 



In the rightmost integral of the right hand side, p(p + E) — p(p) vanishes for E = 0, so we can neglect T in the 
denominator. This gives the final equation: 
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With that approximation, and taking the zero temperature limit, we obtain 



D "" p(E + p) - p(p) dE 

— n F {E) — 

-d-h Po{p) E 



T{T = 0) = (D + p) exp 



1 



JkPo(p) 



dE p (p + E) - p Q (p) 



-D-fj, 



E 



Po(p) 



(A12) 



(A13) 



We note that T(T = 0)/Tk = 2e 7E /it ~ 2.26 for such level of approximation. For finite te mpera ture, we can 
replace n F (E) with 6(—E) in the integral over the density of states in the right-hand side of Eq. (A12|. Using that 
approximation, we can write: 



In 



/T(T = 0) 



du 
4 cosh 2 u 



In 1 



4T 2 u 2 
T 2 



(A14) 



These equations imply that T(T) = T(T = 0)<p(T(T = 0)/T), so given the relation between T(T = 0) and T K , 
T(T)=T(T = 0)<p(T K /T). 



Appendix B: Asymptotic approximation for the Green's function in the presence of warping 

The exact Green's function of the surface electrons is: 

Go(r,iv n ) z 



d 2 k iv„ + p + v F z ■ (a- x k) + | (kl + k 2 _)a z - k . r 



( 2 ^) 2 {iv n +p) 2 -v 2 F k 2 -[^{k\ + k 2 _)\ 



(Bl) 
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Using polar coordinates, we can express the Green's function (Bl ) as a series: 

kdk 



G (r,9,iv n ) 



(2tt) 2 

OO / 

m— — oo \ 



iz + ivpHa+e-^ - as^) + Xk 3 cos{3cj>)a z 



y/( Z 2 + V 2 pk 2 )( Z 2 + V 2 p k 2 + A 2 fc 6 ) 



a ikr cos(0— #) 



2 1.6 



A 2 fc 



\,n\ 



{^(z 2 + v 2 F k 2 ) + ^{z 2 + v 2 F k 2 + A 2 fc 6 )) 2 



( e i64>\m 



(B2) 



where z = v n — ifj,. To obtain an asymptotic expansion of G(r, 9, ii/ n ) to lowest order in A it is enough to consider the 
term m = 0. We then have to consider the integral: 



G (r, 6,iv n ) 



kdk 
(2^ 



iz + ivphjcT+e-* - a-?*) + Afc 3 cos(30)a 2 &ikr cos{(t> _ e) 
^J{z 2 + v\k 2 ){z 2 + v\k 2 + X 2 k 6 ) 



where • • • denotes corrections of higher order in A. After integration with respect to <fi this expression gives: 

„ . „ . f°° 1 AzkJ Q (kr)-v F k 2 (a + e- i6 -a-e i8 )Ji{kr)-i\k 4 cos(?>6)G z Mkr) 

Gn(r, f), iv n ) = / afc , 

io 2^v/(z 2 + v 2 F k 2 ){z 2 + v 2 F k 2 + A 2 fc 6 ) 

s «./ _ ^(a+e-« - ,.eV - £ cos(3^/" . 

Z7T Z7T Z7T 

In the limit A -> 0, we have from Eq. (11.4.44) in Ref. l73l 

I ^ 

Jo z 

k 2 Ji(kr) \v n \-ifi sign(z/ n ) /|i/ n | - i/j, sign(i/ n ) 



kJ (kr) 1 /|i/„| -i,u sign(i/„) 

«fc — ttit = — Ao I r 



Z 2 + Vpk 2 Vp \ <'; 



/' 



(/A- 



Z 2 + TJpfc 2 



'#! 



t'F 



with corrections of order A 2 . We are thus left with the evaluation of the integral: 



I" = 



+ 00 



k A J 3 {kr)dk 



o V(z 2 + {v F k) 2 )(z 2 + (v F k) 2 + A 2 fc 6 ) 



(B3) 

(B4) 
(B5) 

(B6) 

(B7) 

(B8) 



In the expression (B8) we cannot take the limit A — > directly before integrating as we would obtain a divergent 



integral. Instead, we will first make integrations by parts using Eq. (9.1.30) in Ref. [73] to obtain an asymptotic 
expansion. The first integration by parts gives: 



i r +co 

I" = - dkk A JAkr) 

r Jo 



vpk 



vlk + 3A 2 fc 5 



2 r+°o 



dkk 5 J±(kr) 



(z 2 +v 2 F k 2 ) 3 / 2 (z 2 +v 2 F k 2 + X 2 k 6 ) 1 / 2 (z 2 +v 2 F k 2 y/ 2 (z 2 +v 2 pk 2 + X 2 k 6 ) 3 / 2 
1 1 

(Z 2 + V 2 pk 2 f/ 2 {Z 2 + V 2 pk 2 + A 2 fc 6 )!/2 + ( z 2 + u 2 ijfc 2)l/2( z 2 + v | fc2 + ^£6)3/2 



(B9) 



(BIO) 



where • • • denotes the third term in (B9), which will be shown to be negligible in the limit A — ► in App. [C] Using 
another integration by parts, we thus have: 



v 2 r+co 
T" ~ F 



r 2 .10 



k 5 J 5 (kr) 



3v 2 F k 



+ 



v 2 k + 3A 2 fc 5 



2 />+oc 



r 2 .10 



k 5 J 5 (kr, 



(z 2 + v 2 F k 2 ) 5 / 2 (z 2 +v 2 F k 2 + X 2 k 6 y/ 2 (z 2 +v 2 F k 2 ) 3 / 2 {z 2 +v 2 F k 2 + X 2 k e ) 3 / 2 

3(v 2 F k + 3X 2 k 5 ) 



vpk 



(z 2 + v 2 F k 2 fl 2 {z 2 + v 2 F k 2 + X 2 k 6 ) 3 / 2 (z 2 + v 2 F k 2 y/ 2 (z 2 + v 2 F k 2 + X 2 k 6 ) 5 / 2 



(Bll) 



If we take the limit A — > (justification for this can be found in App. O) in the above expression, we obtain a 
convergent integral: 



8vi 



4 r+°° 



k 6 J 5 (kr) _ (\v n \ - in sign(^„)) 3 / \u n \ - i/j, sign(z/„) 



A', 



/ (z 2 + v F k 2 ) 3 vp \ v F 

We can check that the terms proportional to A 2 give contributions that are vanishing in the limit A — > 0. 



(B12) 
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Invoking the asymptotic expressions (B6), (B7), and (B12) together in Eq. (B5), we obtain the following approxi- 



mation for the Green's function in the limit of small warping: 



G(r, 0, iv n ) 



iv-n + M 
2ttv 2 f 



K 



\u n \ - ijx sign (^») 

VF 



r) +i(a + e~ lV -a'e^Kx 



\i/ n \ -in sign( iy) 
vf 



— a A 3 r sign [y n ) cos 30 

v F \ v F J \ v F 



r sign(z/ n ) 



(B13) 



Appendix C: Evaluation of the remainders 



We consider the integral giving the remainder in Eq. ( B9 ) 



q\2 r+°° 

h = — I <n< 



k 9 J 4 (kr) 



r J (z 2 + (vk) 2 ) l / 2 (z 2 + (vk) 2 + X 2 k 6 ) 3 / 2 



With the change of variables k — yup/Au, the integral I3 is rewritten: 

3 r +oc u^hjau) 

3 ~M^rJ Q U {b 2 +u 2 y/ 2 {b 2 +u 2 +u e f/ 2 ' 



(CI) 



(C2) 



with a = r(fi?/A) 1 ' 2 and b 2 = Xz 2 /v F . To obtain an upper bound for the integral, we use 3 successive integration by 
parts to rewrite: 



+00 



u 7 J7(ua)P(u) 



h a 3 Xv F rJ dU (b 2 + u 2 y/ 2 (b 2 + u 2 + u 6 y/ 2 ' 



(C3) 



where the polynomial P{u) has the expression: 

P(u) = 480 u 22 + 1200 b 2 u 20 + 1050 b 4 u 18 ~ 360 u 18 + 315 b 6 u 16 - 3276 b 2 u 16 - 7200 b 4 u 14 - 6084 b 6 u 12 - 432 b 2 u 12 
-1800 b s u 1Q - 432 b 4 u 10 + 1152 b 6 u 8 - 96 b 2 u 8 + 1872 b 8 u 6 - 192 b 4 u 6 + 720 6 10 u 4 + 192 b 8 u 2 + 96 b 10 (C4) 



The numerator in Eq. (C3) is 0(u 57 ' 2 ) while the denominator is 0(u 34 ) for u — > 00 making the integral in (C3) 



convergent. Moreover, an upper bound for the integral is given by: 



+00 



du 



\P(u) 



V2(b 2 + u 2 + u 6 ) 9 / 



(C5) 



,2N-7/2 



< u 7 .We can them majorize the polynomial 



where we have used the inequalities^! | Jt(u)\ < l/\/2 and (b 2 
P(u) by the sum of the absolute value of its monomials. For monomials of degree n > 9, we can also use the inequality 
(b 2 + u 2 + u 6 ) 9 / 2 > u 9 (l + u 4 ) to obtain an upper bound larger than b. For the monomials of degree n < 8, we obtain 
an upper bound of the form: 



+00 



u n du 



{b 2 +u 2 + u 6 ) 9 / 



<6" 



+ OG 



u n du 



Ln+m— 8 



(6 2 +u 2 ) 9 / 2 



+00 



u n du 



(1 + u 2 ) 9 / 2 



(C6) 



From the expression (C4) of P(u), we see that these terms contribute expressions 0(b 2 ). In the case of n = 8, we 
have to consider: 



+00 



u n du 



r-l/b 

• / du 



(b 2 +u 2 +u Q ) 9 / 2 ~ J (1 + u 2 ) 9 / 2 A u(l + u 4 ) 9 / 2 



+00 



du 



(C7) 



Inspecting the polynomial P(u), we see that the contribution of the u 8 terms will be 0(b 2 \ lnfe|). Putting all contri- 
butions together, we see that: 



h< 



< 



a 3 XvFr 



-(C + 0{b 2 \lnb\)) 



Vpr 3 V v F 



(C + 0(A|lnA|) , 



(C8) 
(C9) 
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So 1$ = 0(A 1/2 ) for A — » 0. This establishes that this term gives a subdominant contribution to the Matsubara 
Green's function. 



We can apply the same method to the remainder integrals appearing in Eq. (Bll) 
If we consider the integral: 

X 2 k w J 5 {kr) 



(z 2 +4fc 2 )3/2( z 2 +w 2 nfc2+A 2 fc 6)3/2 



by the same change of variables k = (v/X) 1 ' 2 u, we can transform it into 



h = 



+oo 



u w J§{ua) 



4 v 7 /2 A i/ 2 J (6 2 + U 2 f/ 2 {b 2 + u 2 + u 6)3/2 



Integrating by parts twice, we rewrite: 



1 



du- 



u 7 Jt(uo)P4(u) 



v 7/2 a 2 A l/2 J (&2 + u 2)7/2( fe 2 + u 2 + ^6)7/2 ' 

F 



(CIO) 



(Gil) 



(C12) 



where: 



1,4 ..5 



Pi{u) = 80u u + 100b 2 u l * + 35b* u li + 28u li - 84b z u iL - 204F v;> + 8u y - 92b b u' - 16b 2 u' - 48^ u 



-Wtfu* + 8b°u . 

We can show that the integral remains finite in the limit of A — ► 0, and as a result, 7 4 = 0(A 1/2 ) 
Similarly, for the integral: 



we can rewrite: 



h = 



h = 



+ OG 



X 2 k w J 5 (kr) 



{z 2 + v 2 F k 2 ) 1 / 2 (z 2 + v 2 F k 2 + X 2 k 6 ) 5 / 2 ' 
u 10 J 5 (ua) 



+00 



v V2 A l/2 J (62 + w 2)l/2( 6 2 +u 2 + w 6)5/2 ' 

F 



and by integrating by parts twice, we find: 



1 



+00 



u 7 JT(ua)P^(u) 



2 w 7/ 2 A l/2 J Q (b 2 + U 2 fl 2 {b 2 +U 2 + U 6 f/ 2 ' 



where: 



J.2.,7 



P 5 (u) = 168u i7 + 308 6 2 w i5 + 143 6 4 w lci + 36M lci -108 6 2 u li - 308 6 4 u y + 8u 9 - 164 6 b w 7 - 16 6 2 u 



(C13) 



(C14) 



(C15) 



(C16) 



4„,5 



;6 4 m 



16 6 b u J + 8 6 8 u , 



(C17) 



Repeating the previous reasoning, we again establish that I 5 — 0(A 1/2 ). 

It should be noted that the integrations by part can be iterated as long as the resulting integrals have a finite 
upper bound for small A. This implies that our 0(y/X) estimate for ^3, I4 and i^ is only a conservative one. If the 
integrations by part can be repeated indefinitely, the result of the process is that 13,4,5 = 0(X n ) for any n > 0, a hint 
that the integrals may actually vanish with an essential singularity in the limit of A — > 0. 
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